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We introduce a concept of noncoherent optical pulse discrimination from a coherent (or partially 
coherent) signal of the same energy using a phenomenon of soliton generation. The impact of 
randomisation of the optical signal content on the observable characteristics of solitons generation 
is examined and quantified for a particular example of rectangular pulse. 
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I. INTRODUCTION 



The comparative study of the properties of coherent 
vs noncoherent waves (fields, signals); and ordered vs 
disordered systems plays an important role in both mod- 
ern physics and technological applications (see e.g. [l], Q 
and references therein). The examples range from An- 
derson localization 0] to the interaction of signals with 
noise the latter being one of the major factors limiting 
performance of modern optical communication systems 
(see e.g. [2j]). Optically amplified transmission system 
typically comprise a chain of optical amplifiers (e.g. er- 
bium doped fiber amplifiers) which add noise into the 
system. Amplifier noise accumulated along the fiber links 
degrades the quality of transmission. For instance, sta- 
tistically rare occurrences of parasitic noisy spikes having 
the same energy as the signal pulses might lead to errors 
in the transmission of information. The origin of such 
errors is in applying sub-optimal receivers that make de- 
cision based on information about the energy (or fraction 
of energy) of incoming pulse only and discard any infor- 
mation about the phase (as well as the internal structure 
of the pulse itself). Therefore there is no possible way 
for a standard detection device to distinguish between a 
true logical "one" (characterised by high level of coher- 
ence) and the disordered spurious pulse of the same en- 
ergy. However, the latter "ghost" pulse generally differs 
from the true signal due to a higher level of non-coherence 
and in this Letter we propose a soliton-based device in- 
tended to separate the signal from random fluctuations 
of the same energy. We also discuss here an interesting 
link between the problem of separation of coherent signal 
component from non-coherent noise in fiber optics and a 
theoretical problem of the emergence of localized states 
in a random potential. In particular, we examine the 
generation of coherent localized modes (solitons) from 
random initial field distribution (serving as a potential 
in the corresponding non-Hermitian spectral problem). 

Separation of a signal from noise is one of the funda- 
mental research problems that occurs in a broad range of 
scientific and technological applications. Nonlinear tech- 
niques might be of a special interest in this field as they 
can offer new opportunities that cannot be realized in 
linear systems. For instance, nonlinear optical process- 



ing might enable one to distinguish (in the optical do- 
main) between the signal and noise accumulated within 
the signal bandwidth - something which is not possible 
using linear techniques [^]. Optical fiber nonlinearity is 
of particular interest because the resulting nonlinear sys- 
tems for some range of parameters, can be described by 
well studied soliton models and powerful mathematical 
apparatus can be applied to a variety of physical and 
technical problems (see e.g.0, 0, IE])- Note that optical 
soliton techniques are widely known in the context of 
signal transmission. However their high potential in the 
field of signal processing is yet to be explored. 

Recent progress in advanced modulation formats Q 
and the growing interest in coherent optical communi- 
cations have emphasized the importance of techniques 
that exploit the nature of the signal, in particular, the 
optical phase. Signal processing utilising optical phase 
might offer many advantages over traditional electrical 
processing that relies only on signal intensity. Stable op- 
tical structures such as solitons might play a crucial role 
in the design and development of a new generation of 
optical processing devices that use signal phase as well 
as intensity. Recent demonstration of the quasi-lossless 
fiber span 0] that can be used as one of the building 
blocks of optical soliton devices has rekindled interest in 
soliton-based signal processing. 

In this Letter we propose to use soliton generation as 
a technique to discriminate between coherent (or almost 
coherent) signals and parasitic non-coherent pulses of the 
same energy. We focus here on the fiber-optic applica- 
tions employing the on-off keying data format - when 
logical ones are presented by pulses and logical zeros - by 
empty time slots. We would like to point out, however, 
that the considered approach can be generalised to a va- 
riety of systems and modulation formats. Note that the 
propagation of partially coherent pulses in the Kerr-type 
media has long been a subject of investigation in different 
physical contexts, including random phase modulation of 
temporal solitons [Si , nonlinear Fraunhofer diffraction of 
random fields and other interesting examples [ll[ . 

We would also like to stress that here we are not us- 
ing solitons for transmission of information but rather as 
building blocks of a nonlinear filter element which differs 
the proposed approach from the one employed in con- 
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ventional soliton-based transmission links. The reader 
interested in noise effects in conventional soliton trans- 
mission links is referred to the monograph [3| as well as 
the review 1131. 



II. PROBLEM STATEMENT 

An accumulation of the amplified spontaneous emis- 
sion noise along a fiber-optic link leads to the generation 
of randomly modulated non-coherent background. In the 
on-off keying data format transmission there exist time 
intervals when no deterministic signal is present (time 
slots corresponding to logical zeros). Radiative back- 
ground noise manifests itself as the noise generated ghost 
pulse. It fluctuates and hence can achieve an energy high 
enough to be wrongly detected by a receiver as logical 
one. Note that in the standard approach to signal detec- 
tion, the decision is made taking into account only the 
received pulse energy, and any additional information re- 
lated to the optical phase is not used at all and is lost. 
Any noise field accumulated within the time interval 
and having the optical bandwidth (optical filter band- 
width) B can be presented (following e.g. 0) as a sum 
of statistically independent modes: 



consider here a true soliton system that can be imple- 
mented in lossless fiber spans or in short pieces of 
highly nonlinear fiber with negligible loss effect. Evolu- 
tion in z of the optical field Q{t, z) is then governed by 
the well-known Nonlinear Schrodinger equation (NLSE) 

(Sii). 



dQ 1 d^Q 



dz 



2 dt^ 



(2) 



Here Eq. ^ is written in dimensionless soliton units, 
with time normalized by the pulse width, t t/T, and 



power normalized by Ps = {'^Ld) 



m/ilT^). 7 is 



the nonlinear coefficient in W~^/km, Ld = T^/|/32| is the 
dispersion length measured in km and (32 is the (anoma- 
lous) group velocity dispersion parameter measured in 
ps^/km (see e.g. 0, Q). Energy is measured in units 
Eg — PgT. The particular shape of the input pulse in the 
proposed approach is not critically important, but with- 
out loss of generality, we assume in what follows a real 
rectangular pulse shape with width T < Tb, amplitude 
Qo and energy E — Q^T. 



III. SOLITON GENERATION. 



M/2 
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here int[. . .] means the integer part of an expression. The 
complex coefficients rjk are independent random variables 
with zero mean and variance (rikill,) ~ {2a'^/Tb) Skk' ■ For 
a chain of N optical amplifiers with the same gain coef- 
ficient G; the amplifier inversion factor Usp] and carrier 
frequency of the envelope coq; the noise power spectral 



density per polarization is 2a 



ohujo{G~l) N. The 



average energy of the noise is {Ease) = 2(t^ (M -I- 1), 
which is typically small. However, some rare fluctua- 
tions can produce random spikes with the energy close 
to that of a logical "one" degraded by the same noise. It 
is highly desirable therefore to distinguish between noisy 
(non-coherent) fluctuations and a deterministic localized 
pulse of the same total energy E. Now by normalisa- 
tion of time using a bit length Tb (see below) we can 
express the normalised bandwidth through the value of 
M, therefore, in what follows we will simply quantify the 
bandwidth by the number of modes, M . 

The key idea of the proposed approach is to use soliton 
generation effect for the discrimination between coherent 
and noncoherent pulses. Indeed, the energy of a complex 
Gaussian field ([T|) is characterised by a distribution 
with 2{M+1) degrees of freedom, while the phase of each 
complex harmonic r]k is uniformly distributed between 
and 2tt. Such random phase and amplitude modulation 
typical of parasitic pulses affects soliton generation and, 
thus, can be detected through these changes in a pas- 
sive fiber based device. Without loss of generality, we 



In the case of the real, unmodulated, rectangular, ini- 
tial condition; the energy threshold of the generation of 



A^-th soliton can be readily derived [13|: Ejsf — tt^ (2A^ — 
1)^/4. In the presence of noise, the deterministic formula 



for a number of generated solitons A^^int 0.5 + VE/n 

is replaced by the set of probabilities Pn of generating A^ 
soliton states. Those probabilities Pn{E,(j) will depend 
on the pulse energy, E, as well as the strength of disor- 
der a (and the bandwidth parameter M). Our main goal 
here is to determine how the dependence of an averaged 
number of the generated solitons, (A^) = J2'n=o ^ '-'^ 
input pulse energy evolves when switching from a deter- 
ministic localized input to a noisy non-coherent signal of 
the same energy. 

Here it is pertinent to mention a closely related prob- 
lem of the influence of chirp (both deterministic and ran- 
dom) on soliton generation d, [13, [l^. In particular, 
one may be interested in finding analytical criteria for 
the number of created solitons (and corresponding en- 
ergy thresholds) for different initial pulse shapes. Indeed, 
when the input pulse is a real function waveform a simple 
criterion by Kivshar and Burzlaff (KB) [16] applies that 
relates the number of created solitons with the area of 
the pulse. In a particular case of the rectangular pulse 
one recovers the analytical result above. However, it is 
not straightforward to apply this criterion beyond the 
formal scope of its validity. Indeed, because this crite- 
rion relies only on the area of the pulse it is insensitive to 
any non-trivial phase modulation of the pulse. But any 
such phase modulation (deterministic or random) criti- 
cally affects the soliton content of the pulse - the effect 
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that is missed completely when applying KB criterion. 
For a simple case of a constant linear chirp the prob- 
lem of soliton pulse generation from a Gaussian wave 
form was examined in [l5| where it was shown that lin- 
ear chirp severely decreases the number of soliton states 
in the pulse. It is no surprise then that the KB crite- 
rion is not applicable also to random complex input ([1]). 
As an alternative, a standard WKB approximation can 
be used to estimate the number of created solitons 
but its applicability is restricted to the regime where the 
number of solitons is large. So unfortunately, there are 
no reliable analytical criteria for establishing the number 
of emerging solitons for the output in the form ([T]). 

A. Zakharov-Shabat spectral problem for random 
potentials. 

According to the inverse scattering transform theory 
any initial condition of the NLSE ([2]) evolves into a 
combination of solitons plus quasi-linear oscillating wave 
packets. The parameters of formed solitons, as well as 
their number, can be derived via the non self-adjoint 
Zakharov-Shabat spectral problem (ZSSP) |l7j] which is a 
special case of more general Ablowitz-Newell-Kaup-Segur 
scheme (see [l8|): 

J id%pi/dt + Q%p2 = C^i /o^ 
\ -id^P2/dt + Ri}i =CV'2 ^ ' 

with complex potentials R and Q. The ZSSP spectral 
problem for Eg.© corresponds to the non self-adjoint 
reduction of ([3]): R= ~Q* . Other important reductions 
include R = Q* which corresponds to NLSE with pos- 
itive (normal) group velocity dispersion (or defocusing 
nonlinearity) and the case R = 1 that leads directly to 
the Hermitian Schrodinger equation with the "energy" 
given by connected with the Korteveg-de-Vries equa- 
tion. In any reduction, the random nature of the poten- 
tial leads to interesting properties of the spectra. In par- 
ticular, the reduction which yields the linear Schrodinger 
equation with random potential is known to exhibit An- 
derson localization P, ]^ , which manifests itself as an 
exponential decay of the transmission coefficient through 
a disordered potential with the growth of the width of 
the potential, T. The corresponding decrement t(C) is 
known as localization length (or rather, localization time 
in our notations). For the self-adjoint reduction R — Q*, 
a similar phenomenon exists [l9| where it was proven 
that the localization length does not depend on the spec- 
tral parameter In what follows we will focus on the 
focusing NLSE ^ with non-selfadjoint reduction where 
each (random) discrete eigenavlue located in the upper 
complex half-plane of ^ corresponds to a bright optical 
soliton emerging from a noisy input pulse. As mentioned 
in the previous section the existing deterministic criteria 
for soliton creation are not applicable to the noncoherent 
complex pulses of the form ([1]). Therefore we have to 
resort to the numerical Monte Carlo (MC) simulations of 



the stochastic ZSSP eigenproblem ([3]) and the problem 
of determination of the average number of solitons, (N) 
reduces thus to the calculation of the distribution and 
the average number of discrete eigenvalues of ZSSP. 



IV. CONCEPT OF SOLITON DISCRIMINATOR 
AND RESULTS 

Now we present the results of the MC simulations for 
different system parameters. Those are to be compared 
with the analytical result (as discussed above) for a de- 
terministic rectangular pulse of the same energy. Again 
we would like to convey the idea that we are specially 
interested in rare fluctuations when the energy of the 
"ghost pulse" coincides with that of a real pulse. In the 
MC simulation we had to insure that for each realisation 
of the input field the total energy was fixed at the same 
level, E. This was achieved by picking (M -I- 1) complex 
Fourier harmonics with equal amplitudes (fixed by the to- 
tal energy E) and independent random phases uniformly 
distributed between and 27r. Figure 1 shows the aver- 
age number of emerging solitons for totally non-coherent 
input (H)) versus pulse energy (normalised by the soli- 
ton energy), for different number of random harmonics 
in log-log scale. It is seen from Fig. 1 that the introduc- 
tion of randomness dramatically reduces the probability 
of soliton generation; and the average numbers of solitons 
can be less then unity, even for relatively high values of 
energy. This demonstrates a feasibility of the proposed 
concept of soliton discriminator since it shows that a to- 
tally non-coherent input pulse in unlikely to produce a 
stable soliton content. An interesting question remains 
open of whether the different curves in Fig.l have a fi- 
nite size scaling (ESS) with M, i.e. whether there exists 
an effective "phase transition" in M. In other words, 
whether there exist a universal function / and parame- 
ters a, /3 and Eth so that the average number of solitons 
scales with energy as M"" f{{E - Eth)/M'^). Standard 
algorithms [13] have failed to establish such a scaling, 
but the results are still inconclusive and more advanced 
analysis is needed. 

Another observation is that the deterministic thresh- 
olds of soliton creation (shown in Fig. 1 as vertical dashed 
lines) also do not apply for the disordered input. We de- 
fine a soliton creation threshold for a disordered pulse as 
the first value of E where the average value of solitons 
is non-zero. The computationally resolved thresholds for 
different numbers of harmonics are given by the left end 
points at the graphs in Fig.l. As seen from Fig.l the data 
points in the vicinity of the thresholds tend to fluctu- 
ate which is the natural consequence of the Monte-Carlo 
averaging procedure: the first non-zero contribution to 
the average number of solitons come from the extremely 
rare fluctuations that are notoriously difficult to sample. 
In our simulations we saw no solitons created below the 
first deterministic threshold Ei = 2.46 after 1000 Monte 
Carlo runs. Because ESS algorithm failed to establish the 
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Energy, soliton units 



FIG. 1: (Color online) The average number of emerging soli- 
tons versus pulse energy for non-coherent input pulse for dif- 
ferent number of modes M. Vertical lines show deterministic 
thresholds {Ei = 2.46, E2 = 22.46, E3 = 61.8, and so on) of 
N-soliton generation in soliton units. 
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FIG. 2: (Color online) The probability distribution of the 
number of created solitons (solid red bars) and the Poissonian 
fit (dashed bars). 



unique threshold energy Eth we cannot rigorously prove 
that such processes are forbidden. This issue requires 
a separate study with higher resolution and advanced 
Monte-Carlo techniques (like e.g. multicanonical sam- 
pling) . This however lies beyond the scope of the current 
paper. 

Fig.l shows only the average number of solitons and 
not the corresponding probability distributions of soliton 
formation P/y. It is instructive to study the individual 
probabilities P/v, and in particular verify whether they 
follow a Poissonian distribution. Therefore, in Fig. 2 
we also present separately the probability distribution 
obtained via Monte Carlo simulation of a system with 
M = 8 noisy harmonics and energy E = 360 (in soliton 



units) . The results are compared with the Poissonian fit 
with the same average (N) = 0.73. One can see that 
the numerical distribution is somewhat wider and less 
localized than the Poissonian fit. 

Let us now illustrate different possible scenarios of a 
"ghost pulse" propagation by providing space-time traces 
of pulse evolution for different values of parameters E and 
M. The goal of this detour is to help the reader to gain 
a qualitative understanding of "dangerous" noise fluctu- 
ations that are able to trigger the creation of spurious 
solitons during pulse evolution. In Fig. 3 we provide 
space-time traces of pulse evolution from the stochastic 
input ([1]). We opted here for the following system param- 
eters: Tt = 60 ps, f32 = -20 psVkm, 7 = 3 W'^ km"!. 
The four scenarios shown in Fig. 3 are: one soliton gen- 
eration at high energy and low number of noisy modes 
(Fig. 3(a)), two soliton generation for the same param- 
eters (Fig. 3(b)), no soliton created for high energy and 
high number of modes (Fig. 3(c)) and, finally, no soliton 
generated at sub-threshold low energy (Fig. 3(d)). 

So far we have only considered the evolution and soli- 
ton content of a purely disordered pulse of the form 
dl]). However typically in optical communications a sig- 
nal arriving at the decision device is only partially non- 
coherent and generally can be presented as a determin- 
istic component Qo(^) mixed with an accumulated noise 
q{t), Q{t) = Qo{t) + q{t) where q{t) is given by H]). For 
this case we compute the number of emerging soliton 
states for the ZSSP ^ for the input pulse Q{t) rather 
then for q{t). Figure 4 shows the impact of noise on 
the generation of solitons in this case. The contribution 
of noise to the initial pulse (i.e. the measure of non- 
coherence) is characterised by an optical signal-to-noise 
ratio, OSNR = E/{2a'^ (M -I- 1)) a quantity often used 
in optical communications. Figure 4 demonstrates two 
important limits: large OSNR corresponds to the deter- 
ministic limit described analytically above, and the limit 
of small OSNR leads to the situation of fully non-coherent 
pulse similar to Fig. 1. It also gives a qualitative estimate 
of the level of non-coherent noisy content in the pulse that 
can be distinguished by such a nonlinear soliton-based fil- 
ter. It is seen from Fig. 4 that pulses with OSNR less 
than roughly 10 dB have quite different characteristics in 
terms of the number of generated solitons compared to 
their deterministic counterparts of the same energy (top 
of the picture). This means that if the non-coherent share 
of the total pulse energy content is more than approxi- 
mately one tenth, the proposed device should be able to 
discriminate between such partially coherent signals and 
non-coherent noisy ghost pulses. 

The action of a soliton discriminator is illustrated in 
Fig. 5 where a typical evolution of a given pattern 10100 
is shown in the two cases: (a) when all bits are determin- 
istic and (b) when the central bit is a noisy, non-coherent 
ghost pulse of the same energy as a deterministic one in 
(a). For the particular realization of the noisy input (bot- 
tom panel) we picked M = 64 harmonics with the same 
total energy as the coherent pulse in (a) . Figure 5 demon- 
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FIG. 3: (Color online) Possible scenarios of stochastic soliton creation. The energy (in soliton units) and the number of modes 
are: E = 360, M = 8 for (a) and (h), E ^ 360, M = 64 for (c) and E = 13.5, M = 16 for (d). 
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FIG. 4: (Color online) The averaged number of solitons versus 
the pulse energy E and OSNR. Here M = 64. 
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FIG. 5: (Color online) The evolution of a given pattern with 
the real central "one" (a) and a spurious central "one" (b) 
one of the same energy 



strates how non-coherence prevents the creation of soli- 
ton structures, and one readily observes at the top panel 
(a) of the figure that most of the pulse energy quickly 
settles in an emerging single soliton which is reverse to 
the case (b) at the bottom panel. In the real world units 
this discrimination can be observed at the propagation 
scales of several kilometers for a standard fiber or just 
several meters for highly nonlinear fibers. This implies 
that a potential soliton discriminating device will not not 
require large fiber spans and could thus be quite portable. 



We would like to emphasize that Fig. 5(b) presents just 
a single realization of the event when a noisy background 
accumulates energy comparable to that of the logical one. 
In other words the proposed soliton discriminator might 
be particulary efficient when the system performance is 
limited by such fluctuations. Figure 5 illustrates that the 
proposed device is feasible and is capable of discriminat- 
ing between coherent and non-coherent inputs even if the 
energies are the same and no difference would have been 
detected by a standard integrate-and-dump receiver. 

Let us discuss now an issue of the speed of convergence 
of a disordered input to a final (multi)soliton state. In- 
deed, when analysing a general pulse-to-soliton conver- 
sion problem there are quite a few examples when such 
a convergence to a final sech-shape solution is very slow 
(see e.g. [2l|, [g^, HI]). This, however, does not have a 
serious impact on the proposed concept of soliton dis- 
criminator. Indeed, in order for the soliton discriminator 
to operate one does not have to wait until the ampli- 
tude of the emerged soliton stabilises at a certain fixed 
level. It is not even required at all for the solution to 
stabilise at a prescribed "clean" sech-shape. The impor- 
tant point is to have a solution with a finite energy after 
the radiation has been shed - this will ensure that true 
logical "ones" survive the transmission. The final shape 
of the emerging solution is not of importance for this 
application, it just has to be localized in time domain. 
In the meantime it is desirable that a pulse of the same 
energy, but generated from random Fourier harmonics - 
Eq. ([T]) - will simply collapse and will not give rise to any 
soliton state at all. And Fig. 5 illustrates just that: a 
"ghost one" collapses into linear radiation very quickly 
(and therefore will rightly be detected as "zero") while 
the true logical "one" evolves into an asymptotic soliton 
state. True, the convergence in the latter case is slow, 
but the most important fact is that it survives and still 
carries enough energy to be rightly detected as "one" . 

Evidently, there is a number of technical issues to be 
resolved before any practical implementation of the pro- 
posed soliton discriminating device can be achieved. In 
particular, timing jitter of partially non-coherent signals 
should be accounted for in the analysis of errors. Also 
pattern effects and soliton interaction must be consid- 
ered in any practical implementation. However there is 
an important factor that makes the concept of a soli- 
ton discriminator quite feasible. All the adverse effects 
listed are more pronounced when the propagation dis- 
tance increases (e.g. timing jitter) and cannot be ne- 
glected when propagating through distances in the or- 
der of a couple of dispersion lengths Ld- However, as 
stressed earlier the soliton discriminator can operate at 
much smaller distances. A distance of several kilome- 
ters or even meters usually suffices thus diminishing the 
majority of mentioned detrimental effects. In addition, 
a number of standard well developed control techniques 
(modulation, optical filtering) can be applied to reduce 
these effects even further. 
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V. CONCLUSIONS 

We have proposed soliton generation as a candidate 
technique to discriminate between coherent signals and 
parasitic non-coherent (or partially coherent) pulses of 
the same energy. The principal idea was to use lossless 
transmission spans (modelled by integrable NLSE) to fil- 



ter noisy dispersive part of the solution from the emerging 
soliton part. Using numerical modelling we have demon- 
strated the feasibility of the proposed approach. We have 
also determined how the average number of the generated 
solitons is affected by the presence of noise in the partially 
non-coherent input signal. This work was supported by 
the Royal Society and the EPSRC. 
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